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AN  EXTENSION  OF  SPITZER'S  INTEGRAL 
REPRESENTATION  THEOREM  WITH  AN  APPLICATION 


by 
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D.N.  Shanbhag 

King  Saud  University.  Saudi  Arabia. 
University  of  Pittsburgh.  U.S.A.  and 
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Summary  Using  a  new  approach,  an  extended  version  of  Spitzer's  Integral 
representation  for  stationary  measures  of  a  discrete  branching  process  is  ob¬ 
tained.  This  result  is  used  to  provide  a  complete  solution  to  a  problem  in 
damage  models  satisfying  a  generalized  Rao-Rubin  condition. 


Key  Words;  Branching  process,  Cauchy  equation.  Damage  Models,  Potential  theory 
of  Markoff  chains,  Rao-Rubin  condition. 


AMS  Subject  Classification:  60E05,  62E10. 


1.  INTRODUCTION 


Consider  a  modified  discrete  branching  process  {Zq;  n«0,l,...}  with  one 
step  tansition  probabilities  given  by 


(1.1) 


f  c  p.  •  i*  0,1,. . .  ;j  ■  1,2, . 

PCZ»rt-i1Vl>*  ,  \  a)  ( 

Ll-c+c  pn  ,  1*0,1,...;  j- 


where  0<  c<^l  and  (p^  * }  is  the  i-fold  convolution  of  some  nondegenerate  prob¬ 
ability  distribution  (p^ }  with  itself  for  i>  0  and  the  degenerate  distribution 
at  zero  if  i«  0.  (It  is  seen  that  the  process  reduces  to  Bienaymd-Galton-Watson 
branching  process  when  c*l).  Define 


(1.2)  m-  l  jp  ,  m*  -  £(jlogj)p. 

j-1  3  j-1  3 

In  this  paper,  we  extend  the  Spitzer's  integral  representation  for  stationary 

* 

measures  of  a  Bienaym£-Galton-Watson  process  with  m<  1  and  m  <  °»  to  a  branch¬ 
ing  process  of  the  type  defined  in  (1.1).  Such  an  extension  could  be  arrived  at 
by  using  arguments  similar  to  those  of  Spitzer  (1967)  based  on  the  potential 
theory  of  Markoff  chains.  However,  we  give  a  new  approach  based  on  Bernstein's 
theorem  on  absolutely  monotonic  functions.  We  use  the  extended  version  of 
Spitzer's  theorem  to  obtain  a  complete  identification  of  the  solutions  to  a 
certain  functional  equation  in  damage  models  considered  earlier  by  Talwalker 
(1980)  and  Rao  et  al.  (1980),  and  thus  establish  a  link  between  branching  pro¬ 
cesses  and  damage  models. 


2.  AN  EXTENDED  VERSION  OF  SPITZER'S 
INTEGRAL  REPRESENTATION 


We  establish  the  following  theorems. 

Theorem  1.  The  generating  function  U(s)-  of  any  stationary  measure 

of  the  modified  branching  process  defined  in  (1.1)  is  analytic  for  |s|  <  q 
where  q  is  the  smallest  non-negative  root  of  the  equation  s *  f(s),  and  (if  normal¬ 
ized  so  that  U(Pq)s  1)  satisfies  the  equation 

(2.1)  cU(f (s))  -  c+U(s) 


where  f(s)  is  the  g.f.  (generating  function)  of  {p^}.  Conversely,  if  U(s) 


£n 


s^ ,  I9I  <  d  satisfies  (2.1),  then  (n  }  is  a  stationary  measure. 

13  2 


Proof  If  U(s)  is  analytic  for  |s  <  for  some  Sq >  0,  then 


00  m  00 

U(s)  -  l  n.8^  -In.  I  p.  s 

j-1  J  i-1  1  j-1 


c  l  n.([f(s)]i  -  tf (0) 31) 
i-1  1 

c(U(f(s))  -U(pn)),  |s|<sr 


which  implies  that  (2.1)  is  valid  at  least  for  js|  <  sQ.  Then  using  the  arguments 

exactly  as  in  A.N.  (Theorem  2,  p.  68),  [where  we  use  the  abbreviation  A.N.  for 

the  book  by  Athreya  and  Ney  (1972)],  we  find  that  U(s)  is  analytic  for  |s|  <  q. 

The  converse  assertion  easily  follows  by  equating  the  coefficients  of  s^  in  (2.1) 

* 

Theorem  2.  (An  extended  version  of  Spitzer's  theorem).  If  m<  1  and  m  <  » 
with  f  as  in  Theorem  1,  then  for  every  probability  measure  v  on  [0,1) 


(2.2)  U(s)  -  K  U(s,t)dv(t) 

J[0,l) 

is  Che  g.f.  of  a  stationary  measure,  where 

U(s,t)  -  £  [exp{(8(s)-l)mn  fc}-  exp{-mn  C}]cn  t 

+  * 

with  8(s)  as  the  unique  p.g.f.  among  those  vanishing  at  s  ■  0  and  satisfying  the 
equation 

(2.3)  8(f (s))  -  m8(s) +  1-m. 

Conversely,  every  stationary  measure  has  the  representation  (2.2)  for  some 
probability  measure  v  on  [0,1). 

Proof  The  first  part  of  the  Theorem  is  easily  verified.  To  prove  the 
converse,  it  is  sufficient  to  establish  that  the  representation  is  valid  for  every 
s c  [0,1).  Define  then  for  every  sc  [0,1) 

U*(s)  -  U(8~1(s) ) 

where  8  is  as  mentioned  in  the  statement  of  Theorem  2.  In  view  of  Theorem  1  and 
the  equation  (2.3)  we  have 

(2.4)  c^U  (mns+l-  mn)  *  U  (s)  +  5  ,  n«  1,2,. . . , ;  0  <  s  <  1, 

n  — 

where 

.  n 

(■  c(f^-)  if  c*  1, 

?»-j 

f-  n  if  c-  1. 


We  can  write  (2.4)  also  as 


*For  an  interpretation  of  B  in  branching  processes,  see  A.N.  (p.  17). 
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U*(s)  -  c^U+n^Q  ®_1(s)))- €  ,  n*  1,2, ;  0<s<l, 

n  n  — 

where  (^(s)  “  [fn(3)-l]/mn  with  fQ(s)  as  the  n-th  iterate  of  fCs).  Consequently 

(2.5)  U*(s)  -  lim  {c'hKl+ni^CO)- onQn(0)8n(8_1(s)))-Cn} 

n-*-  “ 

for  0  <  s  <  1,  where  8  (s)  »  [Q  (0)  -  Q_(s)]/Q  (0)  .  If  0<  s.  <  s<  s,<  1,  then  noting 
—  n  nun  —  l  z. 

in  particular  that  8  -*■  8  pointwise  (see  A.N.,  p.  47),  and  8  and  8  1  are  strictly 
increasing  on  [0,1),  we  find  that 

(2.6)  -mnQn(0)8n(8"1(s1)) <  -mnQn(0)s-  -mnQQ(0)8(B“1(s) ) 

<-mnQ  (0)8  (8-1(s,)) 
n  n  z. 


for  large  enough  n.  Since  U  is  increasing  on  [0,1)  we  obtain  from  (2.5)  using 

(2.6) 

(2.7)  U*(s^)  <_  lim  {cnU(l+mnQn(0)-mnQn(0)s) 

n  ■*co 

<_  lim  {cIlU(l+ mIlQn(0)  -  mnQn(0)8)  -  (82)- 

n-*  00 

it  it 

It  is  seen  that  U  (s)  is  continuous  and  U  (0)  ■  0,  and  hence  from  (2.7)  we 
have  for  se  [0,1) 


(2.8) 


U  (s)  «  lim  {c'hjd+m^  (0)  -  mnQ  (0)s)  -£  } 

n  n  n 

n 


-  lim  {cn[U(l+mn 
n  -*■« 


(^(0)-  mnQn(0)s)-  U(l+«n 


<^(0))]} 


for  se  [0,1),  with  the  limits  well  defined.  Since  for  each  n>_l,  the  expression 
within  the  second  limit  in  (2.8)  can  be  expressed  as  the  g.f.  of  some  non-negative 
sequence  for  se  [0,1),  the  extended  continuity  theorem  given  in  Feller  (1971, 
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p.  433)  implies  that  U  is  the  g.f.  of  a  non-negative  sequence.  Define  now  a 
function  U  on  (-®,1)  such  that  its  restriction  to  [0,1)  is  indeed  U*  and  for 
any  se  (-®,0) 

U(s)  -  cnU*(mns+  l-m11)-  K 

n 

where  n  is  the  smallest  integer  for  which  mQs  +  1-  mne  [0,1).  It  is  easily 
seen  that 

(2.9)  cUU(mtls+  l-m11)  -  U(s)+£  ,  n-1,2,...;  se  (-*,1) 

n 

A  ^  A  A 

which,  since  U«  U  on  [0,1),  implies  that  U'  (i.e.,  dU(s)/ds,  se  (-■*», 1))  exists 

A  A 

and  is  absolutely  monotonic  on  (-“,1)  with  lim  U*  (s)  *  0  as  s-*- -®.  Since  U(0)«  0, 
we  conclude  from  Bernstein's  theorem  that  for  some  measure  u  on  (0,®) 

A  f  gx 

(2.10)  U(s)  -  (e  -l)dp(x),  se  (—,1). 

>  (0,®) 

From  the  equation  (2.9),  we  find  that  the  measure  u  is  such  that 

(2.11)  cn  (eSXm  -De^1  m  du(x)  »  I  (eSX-l)du(x) 

J(0,oo)  J(0,®) 

n  *  1,2,...;  se  (-*,1). 


If  Sq  e  (-«,1)  and  se  (-®,0]f  then  subtracting  the  identity  (2.11)  from  the 
corresponding  identity  with  s  replaced  by  s+sn  we  obtain 


(2.12)  c 


_ n  s .  xqi  / .  n. 

sxm  .0  , .  (l-m  ) x  ,  /  . 

e  (e  -l)e  du(x) 


(0,®) 


^(0  ,®) 


eSX(e  0  -l)du(x) 


which  is  valid  for  all  se  (-®,0J.  Then,  in  view  of  the  uniqueness  theorem  for 
Laplace-Stieltj es  transforms,  (2.12)  Implies 


t  srtx  s_xm  ik 

(2.13)  .  .  (e  0  -l)dn(x)  -  c“  ,  (e  0  -l)eU"m  )xd|i<x) 

J [.“.m"'1)  'll,.'1) 

for  n«  0,  +l,+2,....  Since  Sq  in  (2.13)  is  arbitrary,  we  conclude  using  (2.10) 
that 

*  *  00 

(2.14)  U  (s)-U(s)  -  l  ,  (eSX-l)dti(x) 

.  n  n- 1, 
n— [m  ,m  ) 

I  r  n,  -(l-s)xm  -xm  .  x.  .  . 

■  2.  c  (e  -e  )e  du(x) 

J  [l,m  )  n*-«» 

f  “  n-t/,  v  n-t 

-K  {  !  «"■'(.•■  )Mv(t) 

J[0,l)  n— 

for  s£  [0,1),  where  K  is  a  positive  constant  and  v  is  a  probability  measure  on 
[0,1)  such  that  for  every  Borel  subset  A  of  [0,1)  we  have 

(2.15)  v(A).K*1f  .Viog*/10*"  du(x>, 

JSA 

where  S »  (x:  (-log  x  /logm)  e  A}.  [The  operation  of  interchanging  the  order 
of  summation  and  integration  in  (2.14)  is  justified  by  either  Fubini's  theorem 
or  monotone  convergence  theorem.]  The  required  result  now  follows  on  observing 
that  U(s)*U  (8(s)),  s  e  [0,1). 

Remark  1.  In  the  definition  2  on  page  4<32  of  Feller  (19  71),  0,°°  should  be 
changed  to  0,°°  to  make  the  extended  continuity  theorem  given  on  page  433  non- 
ambiguous.  [We  used  this  theorem  in  the  proof  of  our  Theorem  2.] 

Remark  2.  If  we  define  G(i,j)  to  be  the  Green  function  corresponding  to  the 
substochastic  matrix  (p^j)»  i,j*l,2,...,  then  an  argument  essentially  as  in 
A.N.  (p.  70)  implies  that  for  every  subsequence  of  positive  integers  {k  }  such 
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that  the  fractional  part  of  (-log  k^)/logm  +  t  with  te  [0,1)  we  have 


l  G(k  ,  j)s^ 

j-1 


6(^,1) 


-*•  U*(s,t) 


U(s.t) 

U'(0,t) 


where  U(s,t)  is  the  U(s,t)  of  Theorem  2  with  m11  C  replaced  by  (-Q(0)mn  t) ,  Q(0) 
being  as  defined  in  A.N.  (p.  40),  and 


U’(0,t)  -  -aufat-)- 


ds 


s-0 


We  note  that  the  proof  in  AN  (p.  69)  remains  valid  in  the  present  case  as  well 

* 


with  U(s,t)  changed  to  U  (s,t).  Consequently,  the  Poisson-Martin  integral  re¬ 
presentation  for  a  stationary  measure  analogous  to  the  one  given  in  Seneta  (1973, 
p.  151)  for  a  super  regular  vector  yields  the  validity  of  our  Theorem  2  with 
U(s,t)  replaced  by  U(s,t).  Now,  it  is  easy  to  verify  that  Theorem  2  remains  valid 
even  when  Q(0)  is  not  necessarily  equal  to  -1  and  hence  we  have  an  alternative 
proof  of  Theorem  2. 

Remark  3.  The  version  of  Theorem  2  with  U(s,t)  replaced  by  U(s,t)  of  Remark 
2  yields  Spitzer's  theorem  when  c  *  1. 


Corollary,  If  f(s)  ■  1-m+ms,  then  every  sequence  {n^}  is  a  stationary  measure 


iff  it  is  of  the  following  form: 


(2.16) 


I 


J  n“-°o  '  [0,1) 


,  n-t  (n-t)j 
n-t  -m  m  J 


dv(t)  ,  j  -  1,2,.  . . 


where  v  is  a  probability  measure  on  [0,1),  and  K  is  a  positive  constant  as  in 
Theorem  2. 

The  result  (2.16)  is  obvious  from  Theorem  2  in  view  of  the  fact  that 
8(s)  »  s. 


3.  AN  APPLICATION  TO  DAMAGE  MODELS 

Let  (X,Y)  and  (X',Y')  be  2-vector  random  variables  with  non-negative  integer 

valued  components  such  that  X  and  X'  have  the  same  marginal  distribution  {g  } 

n 

with  gg  <  1»  and  Y  and  Y'  are  such  that  for  each  n  with  gn>  0 

P(Y-r|x-n)  -  C“)TTr(l-ir)n“r,  r=0,l,...,n, 

P(Y'  -  r'  |X'  -  n)  -  (”)tt  'r(l— it  ’  )n  r,  r*0,l,...,n, 

where  0<  it',tt<  1  are  fixed.  Talwalker  (1980)  and  Rao  et  al.  (1980)  considered 
the  problem  of  characterizing  the  distributions  of  X  and  X'  by  the  equation 

(3.1)  P(Y»  r)  -  P(Y'- r|x'«  Y'),  r=0,l,... 

which  is  an  extended  version  of  the  Rao-Rubin  condition  (Rao  and  Rubin,  1964) . 

It  is  seen  that  (3.1)  is  equivalent  to 

(3.2)  G(1-tt  +  tts)  •  G(stt*)/G(it’)  ,  | s  J  <_  1 
where  G  is  the  g.f.  of  {gnh 

When  tt=*  tt',  (3.1)  is  the  Rao-Rubin  condition,  and  it  is  shown  by  Rao  and  Rubin 
(1964)  using  Bernstein's  theorem  and  in  a  simpler  way  by  Shanbhag  (1974),  that 
in  this  case  G  is  the  p.g.f.  of  a  Poisson  distribution.  When  tt  >  tt',  it  is  easily 
seen  that  G  is  the  p.g.f.  of  a  binominal  distribution  with  an  arbitrary  index  and 
success  probability  (tt— tt  ' )  / tt(1— it  '  )  . 

When  tt  <  tt',  the  picture  is  totally  different,  and  the  family  of  distributions 
{gn}  for  which  (.3.1)  holds  is  somewhat  curious  and  fairly  large.  Talwalkar  (1980) 
and  Rao  et  al  (1980)  identified  the  family  as  a  mixture  of  Poisson  distributions 
with  the  mixing  measure  itself  satisfying  a  further  functional  equation.  The 


following  theorem  gives  a  complete  identification  of  the  solution  and  provides 
a  more  satisfactory  answer  and  a  rigorous  proof  to  the  characterization  problem 
than  that  given  earlier. 

Theorem  3.  Let  (X,Y)  and  (X',Y')  be  2-vector  random  variables  as  considered 
above  with  0<  ir<  ir*<  1.  Then  (3.1)  is  valid  if  and  only  if 


"  K 


l 

—CO  . 


\  n-t  -(ir/ir*)11  C  (tt/tt * ) 


[0,1) 


^,(l-TT)^dV(t)  ’  J  "°  ’1»- 


where  v  is  a  probability  measure  on  [0,1)  and  c  is  a  real  number  lying  in  (0,1) 
and  K  is  a  normalizing  constant. 

The  result  follows  from  the  Corollary  to  Theorem  2  (see  equation  (2.16)) 
since  (3.1)  is  equivalent  to  (3.2)  which  can  be  written  as 


with 


cU(l-m+ms)  =  U(s)+cU(l-m) 


io’irjW1-").  »-? 


»<•>  -  s>  -  «o- 


Observe  that  (3.2)  implies  that  the  p.g.f.  G(s)  is  defined  for  all  s  such  that 
|s|  <  it'  (1— tt)  / (-TT *  — tt)  and  it  satisfies  (3.2)  for  every 

f  TT  *  (1— TT)  TT'd-Tr)^ 

l-  /  •  \  ,  »  /• 

V7T  —  n  )  TT  — TT 

Remark  4.  If  the  measure  v  in  Theorem  3  is  taken  as  the  Lebesgue  measure 
on  [0,1),  then  the  distribution  {g  }  in  question  reduces  to  a  negative  binomial 


distribution. 


.0 


Remark  5.  It  is  interesting  to  note  that  if  a  p.g.f.  G  satisfies  (3.2) 
simultaneously  for  two  pairs  (it^tt*)  ,  i-  1,2,  where  0<  iri<  tt^<  1  and  (log  u^-log  tt^) / 
(log  ^“log  irp  is  irrational,  then  G  is  the  p.g.f.  of  a  negative  binomial  distri¬ 
bution  of  the  form 

G(S)  =  {  [ir £(1— ir1)/(-rr;J>-ir1)  ]-a>« 

for  some  a>  0.  Since  the  condition  implies  G  to  be  well  defined  also  on  (1,Sq) 
where  s^3*  ir|(l-ir^)/(irj-ir^) ,  the  result  in  question  follows  as  a  corollary  to  the 
result  of  Marsaglia  and  Tubilla  (1975)  by  noting  in  particular  that  f(x)  « 

G(SQ-Sge  X)/G(0),  x^O  is  well  defined  and  satisfies  the  equation  f(t^+x)  = 
f(ti)f(x),  x  >_  0,  i“  1,2  with  t^  =  log(irJ/ir  ) ,  i"  1,2.  The  same  result  was  es¬ 
tablished  in  Rao  et  al  (.1980)  by  a  different  and  slightly  more  involved  method. 
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